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Abstract 

Newton’s inequalities cjj > c n -ic n +i are shown to hold for the normalized coefficients c n of the 
characteristic polynomial of any M- or inverse M-matrix. They are derived by establishing first an 
auxiliary set of inequalities also valid for both of these classes. They are also used to derive some new 
necessary conditions on the eigenvalues of nonnegative matrices. 


1 Introduction 

The goal of the paper is to prove a conjecture made in : 1] about a set of inequalities satisfied by (the 
elementary symmetric functions of) the eigenvalues of any M- or inverse M-matrix. 

Let (n) denote the collection of all increasing sequences with elements from the set {1,2,..., n.}, let 
denote the size of the sequence a, and let a' denote the complementary or ‘dual’ sequence whose elements 
are all the integers from {1,2,... ,n} not in a. Given a matrix A G <0 nxn , the notation A(a) (A[a]) will be 
used for the principal submatrix (minor) of A whose rows and columns are indexed by a. By convention, 
A[0]:= 1. 

A matrix A is called a P-matrix if A[a\ > 0 for all a £ (n). A is called a (nonsingular) M-matrix if it 
is a P-matrix and its off-diagonal entries are nonpositive. If in this definition the positivity of all principal 
minors is relaxed to nonnegativity, one obtains the class of all M-matrices, including the singular ones. The 
class of inverse M-matrices consists of matrices whose inverses are M-matrices. The M-matrices are an 
important class arising in many contexts (see, for example, (2J Chapter 6]). 

Given a matrix A , let Cj (A) denote the normalized coefficients of its characteristic polynomial: 



The inequalities 

c 2 j (A)>c j - 1 {A)c j+1 {A), j = 1,... ,n — 1, (1) 

are known for real diagonal matrices, i.e., simply for sequences of real numbers (see m and references 
therein), as was first proved by Newton. Since the numbers c, are invariant under similarity, Newton’s 
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inequalities 0 also hold for all diagonalizable matrices with real spectrum, and therefore also for the closure 
of this set, viz. for all matrices with real spectrum. 

It was conjectured in |1’ that Newton’s inequalities are also satisfied by M- and inverse M-matrices 
(and by matrices similar to those). The next section contains proofs of several results on M-matrices and 
symmetric functions culminating in the proof of this fact. 


2 Proof of Newton’s inequalities 


Let us begin by establishing a set of auxiliary inequalities first. Given an n x n-matrix A and nonnegative 
integers mi, m 2 , k, define functions S miiTO2! t as follows 

S mi ,m 2 ,k(A):= J2 A[a]A[p}. ( 2 ) 

cxG{n) ,#a = m 1 , 

0 G (n ) , # P= m 2 , # a fl (3 = k 

Theorem 1 For any M- or inverse M-matrix A of order n and nonnegative integers m < n, k < m, 

(A)/S m 

,m,k {In) ^ l,k (A)/5 ro+1 , m _ llfc (7 n ), (3) 

where I n denotes the identity matrix of order n. 


Proof, by induction. 

Case 1 (induction base). If k = 0, n = 2m, then 0 is a special case of Theorem 1.3 from ?5J. Indeed, 
since n = 2m, the functions and S ro +i, m -i,o are immanants, A :=(m,m) and /u:=(m + l,m — 1) 

are partitions of n, and /i majorizes A. Then the normalized immanant corresponding to fi does not exceed 
the one corresponding to A (beware a typo in 0 where the sign is reversed). If an M-matrix A is non¬ 
singular, then j4 -1 [a] = A[a !]/detH (see, e.g., J3] Section 1.4]), hence S rri!m fi{A~ 1 ) = S m , m fl(A)/(det A) 2 , 
<Sm+i,m-i,o(- | 4 -1 ) = <S'm+i,m-i,o(^4)/(det A) 2 , so the inequality 0 holds for the matrix A -1 as well. 

Now assume 0 holds for all M- and inverse M-matrices of order smaller than n. 

Case 2 (induction step of the first kind). Suppose 2m — k < n and A is an M- or inverse M-matrix. 
Then both normalized functions S m ,m,k{A)/S m ,m,k{In) and S , m+lim _ life (A)/ S m +i,m-i,k{In) can be obtained 
by first averaging the terms AL[a]A[/3] over submatrices of order n — 1 and then taking the average of the 
obtained n quantities: 


S m ,m,k (A) 
Sm,m,k (In) 


*$Vn+1 ,m—l,k {A) 
Sm-\-l,m—l,k (In) 


1 E 

qG (n) ,#at=n—l 

1 E 

n 1 ' 

qG (n) ,#ot=n—l 


Srn ,m,k {A(a)) 
$m,m,k {In— l) 

1 ,m — l,k {A(a)) 

^m+1 ,m—l,k (In- 1) 


But principal submatrices of M- (inverse M-) matrices are again M- (inverse M-) matrices (0 p.113, p. 119]), 
therefore the inductive assumption holds for all submatrices A(a), ffa = n — 1. This implies 0 for the 
matrix A itself. 

Case 3 (induction step of the second kind). Let 2m — k = n and k > 0. First assume A is a nonsingular 
M- or inverse M-matrix. Switch to the dual case: Each Al[a]AL[/3] in the right-hand side of 0 equals 
A~ 1 [a']A~ 1 [j3 l ]/(det A) 2 , the index sets a' and (3' do not intersect, and ffa' + #/3' = 2 (n — m) < n. Hence 


Sr, 


<(A) = 


Sn—m,n—m,o(A ) 

(det A) 2 ’ 


l,k 


(A) = 


m+l,n— m— 1,0 (A- 1 ) 

(det A) 2 


and the functions S'„_ m>r i_ m> o(^4 _1 ), S , n _ m +i, n _ m _i ! o(Al _1 ) are as in Case 2 above. Thus 0 holds for the 
matrix A~ x and hence for the matrix A. So, the induction step of this kind is now proved for nonsingular M- 
matrices and their inverses. But the set of all M- matrices is the closure of the set of nonsingular M-matrices 
(see, e.g., 0 p.H9]), which justifies the induction step for singular M-matrices as well. 
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With all possible cases considered, the theorem is proved. □ 

Now, the theorem can be used to replace Newton’s inequalities by a stronger (but simpler) set of quadratic 
inequalities in the variables A [a]. 

Lemma 2 Let m € {1,..., n} be fixed and let t(m) be the column vector 

t(77l) . {t a ) a£{n) ,#a=m' 

Let denote the Hermitian form 

m 

t(m) i—> t(m)*^ m t{m) := (m(n — m) — (m + l)(n — m + 1)———-—) t a ti 3 . (4) 

U rn-j + l #a ^ =m 

#«n f>=j 

If'km, is nonnegative definite, then the mth Newton’s inequality m holds. 

Proof. Expanding both sides of the mth Newton’s inequality yields 

cl(A) 

Cm —1 1 (-^) 

So, the mth Newton’s inequality is equivalent to 


m / \ 2 

= y^Srn,mj04)/( J i 

j=0 V J 

= § Sm+1 ’ m “ 1,j ( A) /(m + l)(m-l) 


m 


( n-m)'^2 Sm,m,j(A) > (to + l)(n - to + 1) ^ S m +i,m-i,j{A). 
j=o i =0 


(5) 


On the other hand, straightforward counting gives 

(In') 


n \ ( n — J \ ( n — TO 
j) \m - j) \m - j 
n\ / n — j \ fn — m + 1 
j) \m - j - l) \m - j + 1 


hence the inequalities 0 are equivalent to 

(jfl > (to j 1 )*Sm+l,m—lj(^d)- 

Thus, upon replacing each SVn+i,m-ij in the right-hand side of 0 by S m ,m,ji one obtains a set of 

inequalities stronger than Newton’s. Precisely, these stronger inequalities assert that 


m-j 


y^(m(n - to) - (to. + 1 )(n - to. + 1)- j^—^)S m . rn .j > 0, 


j=0 


TO — J + 1 


or, recalling the definitions of SVn,mp and of ’Em, 

a(m)* m a(rn) > 0 where a(m) :=(A[a]) ae ( n } >#CK=TO . 

So, if 'I'm is nonnegative definite, then the mth Newton’s inequality is satisfied. 
Thus, it remains to prove the following. 


□ 


Lemma 3 With the notation of Lemma [ U t{m)*'^ m t{m) > 0 for all t(m) and all m = 1, ..., n — 1. 
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Proof. Consider first the Hermitian form 


5> m : t(m) t(m)*$ rn t(m) := V'j 


t a tp . 


i=o 


# a = #f3=m 

#«np=j 


The representation matrix 


of this Hermitian form is the Gramian, with respect to the standard inner product, for the system of vectors 
( v a ) a where 


I 1 if i £ a 
\ 0 otherwise, 


hence is nonnegative definite. Moreover, the vector e of all ones (of appropriate length) is an eigenvector of 
<f> m . Now consider a form 


: t(m) 


t(m)*<&mt(m) :=y^(m - j + 1) 
3=0 


#c = #/3=m 

#anp=j 


Its representation matrix is obtained by subtracting <E> m from a positive multiple of the Hermitian rank- 
one matrix ee* (precisely (to + l)ee*), therefore all eigenvalues of 4> m are nonpositive except for the one 
corresponding to the eigenvector e, which is strictly positive. Therefore, by [I], the Hadamard inverse 
of the matrix <f> m , i.e., the matrix 

(_i_ 

\m — ffa n/3 + 1 

is nonnegative definite. Finally, \1 > m is obtained from (m + l)(n — m + 1)41 m by subtracting the rank-one 
matrix ee* this time multiplied by (n + 1). The eigenvalue of corresponding to e is equal to zero, since 



e 'i’me = min - m) S m ,mAin) — (m + l)(n — m + 1) N- : -- 

j =0 3=0 J 


Sm,m,j (In) 


= m(n — m)^2 S m , mj (I„) - (m+ l)(n - m + 1) ^ S m +i,m-i,j(In) = 0. 
3=0 j =0 


All the other eigenvalues of I'm are nonnegative, so is nonnegative definite. 
This lemma finishes the proof of Newton’s inequalities. 


□ 


Theorem 4 Let A be similar to an M- or inverse M-matrix. Then the normalized coefficients of its char¬ 
acteristic polynomial satisfy Newton’s inequalities W- 

Also note that a by-product of Lemma [3] is a binomial identity: 

Corollary 5 Ej=oM n ~ m) - (to + 1 ){n - m + l) J!~j 1 )(7) (7-7) = °- 


3 Newton’s inequalities and the inverse eigenvalue problem for 
nonnegative matrices 

As possible applications of Theorem@]one can envision eigenvalue localization for M- and inverse M- matrices 
as well as inverse eigenvalue problems. In the rest of the paper the focus will be on the latter problem for 
nonnegative matrices. 
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The nonnegative inverse eigenvalue problem (NIEP) is that of determining necessary and sufficient con¬ 
ditions in order that a given n-tuple be the spectrum of an entrywise nonnegative n x n matrix. For details 
and history of the problem, see 0 , 0 , Ha, and references therein. 

Two known necessary conditions that an n-tuple A :=(Ai,..., X n ) be realizable as a spectrum of a non¬ 
negative matrices are formulated in terms of its moments 

n 

»»(A):=£A$, 
i=i 

viz. 

Sfc > 0, all k (6) 

sf <n rn ~ 1 s km , all k, to. (7) 

The condition m follows simply from the fact that tr (A k ) is the kth moment of the eigenvalue sequence of 
A, while the condition 0 is due to Loewy and London DU and, independently, Johnson 0. 

Newton’s inequalities proven above result in a third set of conditions necessary for realizability of a given 
n-tuple. Precisely, if A = (Ai,..., X n ) is the spectrum of a nonnegative matrix A and Ai = max |A| is its 
spectral radius, then the set (0, Ai — A 2 ,..., Ai — A„) is the spectrum of an M- matrix X\I — A and should 
therefore satisfy Newton’s inequalities 0 . 

Newton’s inequalities are independent of m and 0 . First of all, it is clear that 0 and are indepen¬ 
dent: for example, the triple (1, —1, —1) does not satisfy 0 but its shifted counterpart (0, 2,2) satisfies 0 , 
while the triple (x/2, i, —i) satisfies 0 but the corresponding shifted triple (0, a/2 — i,y/ 2 + i) does not 
satisfy 0 . 

Moreover, neither the two conditions 0 and 0 together imply 0 nor the two conditions 0 and 0 
together imply 0 . 

Indeed, the conditions 0 and 0 can be satisfied while the conditions 0 may fail. To show this, consider 
the 10-tuple A :=(3,1,1,1,1,1, —2, —2, —2, —2). Its first and third moment are equal to zero, while the rest 
are positive. Now, let us introduce its perturbed version A t :=(3 + ti, 1 + t 2 , 1,1,1,1, —2 + f 3 , —2, —2, —2), 
where the i’s are real and 


t\ + t 2 + £3 >0 

(3 + fi) 3 + (1 + t 2 f + (-2 + f 3 ) 3 = 20, 

which is always possible according to the Linearization Lemma nn p.163], since the system 

t\ + t 2 + t;i >0 

9ti + t 2 + 4t 3 = 0 

is solvable arbitrarily close to the point (0,0,0). The first moment of A t is thus positive, while the third is 
still zero. All the other moments remain positive if (ii, t 2 ,ts) is sufficiently small. So, 0 is satisfied. The 
Newton conditions 0 are satisfied as well, since A t is real. But the condition 0 with k = 1, m = 3 fails. 

To construct an example where 0 and 0 are satisfied but 0 fails, consider the sequence of zeros of the 
polynomial p(x) = x 6 — 6a: 5 + 14a: 4 — 20 .t 3 . It does not satisfy 0: This polynomial is obtained by cutting 
the expansion of (x — l) 6 , whose coefficients satisfy 0 with strict equalities, and then decreasing slightly 
the value (originally 15) of one coefficient. Then the second Newton’s inequality fails. The non-zero roots of 
p are approximately 3.6702 and 1.1649 ± 2.0229L By shifting back by the largest absolute value ss 3.6702, 
one obtains the 6-tuple A :=(a, a, a, 0, 6, b) with a ss 3.6702, b ss 2.5054 + 2.0229L It is not hard, though a bit 
tedious, to check that A satisfies 0 . Since si(A) > 0, this also implies that all moments of A are positive. 
This shows that o cannot be derived from © and 0 . 

In the case the first moment of an n-tuple is zero, Laffey and Meehan [0 established another necessary 
condition, viz., 

(n — l)s 4 > s\. 
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It is also not implied by ®, 0 and 0. An example is provided by the 5-tuple (3, 3, —2, —2, —2). 

Note, however, that the condition 0 with k = 1, m = 2 is exactly equivalent to the Newton inequality 0 
for j = 1. 
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